Apparent Clustering of Intermediate-redshift Galaxies as a Probe of Dark Energy 
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We show the apparent redshift-space clustering of galaxies in redshift range of 0.2-0.4 provides 
surprisingly useful constraints on dark energy component in the universe, because of the right balance 
between the density of objects and the survey depth. We apply Fisher matrix analysis to the the 
Luminous Red Galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS), as a concrete example. 
Possible degeneracies in the evolution of the equation of state (EOS) and the other cosmological 
parameters are clarified. 
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The dark energy, such as the cosmological constant, 
has now turned out to be a necessary element to under- 
stand our universe. There are many indirect suggestions 
for the dark energy, including the age of the universe, the 
formation of the large-scale structure, the number count 
of the galaxies, and so on |Q. More striking evidences 
stem from the combination of acoustic peaks of temper- 
ature fluctuations in cosmic microwave background radi- 
ation (CMB) |g] and the Hubble diagram of the type la 
supernova q|. Natural expectation for the dark energy 
is arisen from the vacuum fluctuations of quantum fields, 
although the smallness of its observed value is extremely 
unnatural [Q. 

Since the energy density of the cosmological constant 
is constant in time by definition, one needs an extremely 
suspicious fine-tuning of 120 digits to set a correct value 
at an initial stage (probably Planck time) so that the en- 
ergy density of the cosmological constant should be com- 
parable with the matter density today. Alternatively, 
such stringent fine-tuning should be moderated if the 
dark energy component consists of some dynamical field 
such as "quintessence" In general, such dynamical 

dark energy affects cosmological observations through a 
time-dependent EOS, which is characterized by a param- 
eter w(t) = p/ p. If w is found to be different from —1, 
the dark energy is proven to be different from the simple 
cosmological constant. Current observational data are 
consistent with the cosmological constant |) , although 
they are still not enough to impose sufficient constraints 
on the EOS. 

Since the dark energy has profound implications on the 
nature of the universe, it is of great importance to explore 
the origin of dark energy Near-future observations will 
enable us to clarify whether the dark energy is a mere 
cosmological constant or something else. Cosmological 
observations, such as the Type la supernovae J|, CMB 
fluctuations ||, |l0), cluster mass function [[ll], weak 
lensing field 113] , and so on, are useful tools to probe the 
nature of the dark energy. 



It is shown [[13| that an application of the Alcock- 
Paczynski (AP) test (lj] to the redshift-space correla- 
tion function of the high-redshift objects can be a useful 
probe of the cosmological constant. In literatures, this 
method is mainly applied to the Lyman-a forest |l5| , the 
Lyman-break galaxies ]l6[j , and the quasars |17|[ , etc. 

A drawback of the high-redshift objects is their sparse- 
ness. Relatively high shot noise prevents to accurately 
determine the correlation function. On the other hand, 
a typical sample of normal galaxies is dense enough to 
have sufficiently small shot noise, while it is too shallow 
to apply the AP test. Instead, a sample of intermediate- 



galaxies, such as the LRGs 18 in the SDSS has a right 
balance between density and depth |l^, |2(| . 

Therefore, a natural expectation is that the cluster- 
ing of intermediate-galaxies can distinguish various dark 
energy models quite well. In fact, the AP test around 
redshift 0.5 is suggested to be promising pl| . The pur- 
pose of this Letter is to quantitatively investigate what 
constraints can be imposed on the parameter space of 
dark energy models from intermediate-redshift galaxies 
and to clarify possible degeneracies with other cosmolog- 
ical parameters. 

The clustering of objects in redshift space is charac- 
terized by the correlation function. In linear regime, the 
full information on the clustering is contained in the two- 
point correlation function, if the non-Gaussianity of the 
density field is negligible. In a homogeneous, isotropic 
space, the two-point correlation function £(r) is a func- 
tion of only a separation r of the two-points. However, 
the redshift-space clustering in reality is neither homoge- 
neous nor isotropic. Peculiar velocities, evolutionary ef- 
fects, and cosmological geometry introduce inhomogenc- 
ity and anisotropy in the observed space. Therefore, the 
two-point correlation function in apparent redshift space 
is generally expressed as £(zi, Z2,9), where z\, Z2 are the 
redshifts of the two-points and 9 is the apparent angle 
between two-points from the observer. The full analytic 
expression for the apparent correlation function in linear 
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regime is known |22| . 

It is not trivial how one can analyze the three dimen- 
sional function £(zi, Z2,&). Instead of directly evaluate 
the two-point correlation function, we have been develop- 
ing a likelihood method for the density field [|l9| pi[ , 
which directly analyze the probability of the observed 
density field given a cosmological model. A fast algo- 
rithm to compute the correlation matrix is crucial in this 
analysis. A detailed description of one of our algorithms 
is given in pej , in which distant-observer approximation 
is assumed. We generalize this algorithm to the one with- 
out this approximation, using the analytic expression of 
p2j . Apart from this generalization, and from inclusion 
of the dark energy degrees of freedom, the algorithm we 
employ in this Letter is similar to that in [^o| , the details 
of which will be published elsewhere. 

Here we briefly summarize major points of the gener- 
alization. In our likelihood analysis of the redshift-space 
clustering, first we place smoothing cells in redshift space 
and then count the number of object rii in each cell i. 
When the non-Gaussianity of the distribution is negligi- 
ble, the correlation matrix dj = (niiij)/((ni)(nj)) — 1 
fully characterizes the statistical property of the cluster- 
ing. The correlation matrix is given by numerical inte- 
grations of the analytic form of Z2, 0) which is gen- 
eralized to take into account nontrivial the EOS for dark 
energy component. Both cosmological distortions and 
peculiar velocity distortions are taken into account in 
this analytic form. A spherical cell approximation pc| ] 
greatly simplifies this integration without implementing 
the numerical integration cell by cell. Then shot noise 
term is added to the correlation matrix. 

Most of the modification of the analytic form of the 
correlation function is originated from the non-standard 
time-dependence of the Hubble parameter. When the 
dark energy has a non-trivial EOS as a function of red- 
shift, w(z), the time-dependent Hubble parameter is 
given by 



H(z) = H Q 



(l + z) 3 fi M0 -(l + ^) 2 ^ 



Kl] 



cxp 3 



I + w(z) 



1 



dz n 



'QO 



1/2 



(1) 



where Qmo> ^Qo are the present density parameter of 
matter and dark energy components, respectively, and 
Oko = ^mo + ^qo — 1 is the present curvature param- 
eter. Evaluation of the analytic form of the apparent 
correlation function requires the comoving distance x( z )i 
the growth factor D(z), and the logarithmic derivative 
of the growth factor f(z) = d\nD/d\na as functions of 
redshift. All these quantities are determined by the time- 
dependent Hubble parameter. The comoving distance is 
simply given by x( z ) = Jo H~ 1 (z)dz. The growth fac- 
tor is the growing solution of the differential equation, 
D + 2HD - §# 2 ^mo(1 + zfD = 0, where dot repre- 



sents the differentiation with respect to the time t. It is 
useful to rewrite this equation into the following set of 
equations, 
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where a = (1 + z) 1 is the scale factor of the uni- 
verse-, and n M (z) = H 2 (l + z) 3 Q M o/H 2 (z), Qq(z) = 
Hq 2 exp(3 Jp j+^dz)flQo/H 2 (z) are the time-dependent 
density parameters of matter and dark energy, respec- 
tively. The Runge-Kutta integration of the set of equa- 
tions (H) , (||) simultaneously gives the growth factor and 
logarithmic derivative of the growth factor. Equations 
(0), (||) are valid even when w evolves with time. 

Once the correlation matrix is evaluated, it is straight- 
forward to obtain expected bounds one can impose on 
a set of model parameters by a given data set, thanks 
to the Fisher information matrix. In linear regime, the 
distribution of the counts n,; is multi-variate Gaussian. 
In this case, the Fisher information matrix has a simple 
form F aP = ^{C^C^C^C.n) @ ||, |§, where C 
is a theoretical model of the correlation matrix, a, (3 are 
indices to indicate a kind of model parameters and C. a 
is a derivative of the correlation matrix with respect to 
a model parameter. An inverse of a Fisher matrix gives 
an estimate of the minimum error variance of a given set 
of model parameters |2f| . 

We consider the LRG sample of the SDSS as a specific 
example. The survey area is assumed to be pi steradian, 
the redshift range to be 0.2-0.4. The number density of 
the LRGs is approximately homogeneous and is given by 
10 -4 /i 3 Mpc~ 3 To reduce the computational cost, 



we set a sub-region of a 10 



314 square degree field 



in a redshift range 0.2-0.4, and fill spherical cells of ra- 
dius 15ft- _1 Mpc in this region. We assume nonlinear ef- 
fects of velocity distortions are erased by this choice of 
smoothing radius. In actual applications, it is important 
to test the choice of smoothing radius to ensure non- 
linear distortions do not bias the result. We use the 
cubic closed-packed structure which has the maximum 
spatial filling factor of 0.74 that can be filled with spher- 
ical cells without overlapping each other (there are other 
ways to achieve this maximum factor such as the hexag- 
onal closed-packed structure, etc.). As a result, about 
2500 cells are placed in this sub-region. The Fisher in- 
formation matrix is scaled according to the ratio of the 
volume of the sub-region to the total volume. 

The model parameters for the dark energy are Oqo and 
w(z). To characterize an evolutionary effect of EOS, we 
employ a parameterization w(z) = wo + WiZ, which is a 
good approximation when the evolution of EOS is mild 
in our interested redshift range of z < 0.4. Therefore, 
we have three model parameters 51qo, Wo, and w\ for the 
dark energy component. There are other cosmological 
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parameters to compute the model correlation matrix. We 
assume a cold dark matter (CDM) power spectrum which 
depends on the shape parameter T = f^Mo^j an d the nor- 
malization erg of mass. The apparent redshift-space corre- 
lation function depends on flMo and bias parameter 6, as 
well as dark-energy parameters and the power spectrum. 
In the following we choose the curvature parameter f2xo 
as an independent parameter instead of Omo- Thus we 
have seven model parameters. In the following analysis, 
we set fiducial values J7qo = 0.7, wq = —1 (or —0.5), 
w'i = 0, f^Ko = 0, h = 0.7, erg = 1, b = 2, and see how 
one can constrain the parameters around this models. 

First, we consider a situation in which only two pa- 
rameters of the dark energy are constrained, assuming 
all the other cosmological parameters are fixed. Fig. 
shows the expected error bounds on two of dark energy 
parameters are varied. In Qqq-wq and wq-wi panels, two 
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FIG. 1: Expected error bounds of dark-energy parameters. 
One of three parameters are fixed. Inner ellipses represent 
the 1<7 uncertainty level of one-parameter probability distri- 
bution. Outer ellipses represent la of the joint probability 
distribution. Solid lines: bias parameter is fixed. Dotted lines: 
bias parameter is marginalized over. 

cases with fiducial values wo = — 1, —0.5 are shown, and 
in flQQ-wi panel, only wq = — 1 case is presented. In- 
ner and outer ellipses correspond to expected ltr bounds 
of one-parameter distributions, and la bounds of two- 
parameter joint distributions, respectively. Solid ellipses 
represent the case in which bias parameter is fixed, as 
well as other cosmological parameters. Since the bias pa- 
rameter is the most important uncertainty in analyses of 
the galaxy clustering, we also plot the expected bounds 
marginalized over bias parameter (dotted ellipses). This 
corresponds to the case when bias should be simultane- 



ously constrained together with dark-energy parameters 
in a same data set. 

One may suspect that the constraints on dark energy 
are imposed mainly from the evolution of the growth fac- 
tor D(z), rather than the AP effect. In which case, the 
uncertainty in a possible bias evolution could demolish 
the constraints. To see this is not the case, we run a mod- 
ified code to calculate the Fisher matrix for an imaginary 
model in which the growth factor does not evolve at all. 
The resulting values of Fisher matrix are only 4%, 30%, 
and 35% smaller than the original values for f^Q-f^Q, wo- 
rn®, and w\-wi components, respectively. Therefore, the 
AP effect does mainly constrain the dark energy param- 
eters. 

In the f2Qo-u>o panel, non-evolved EOS models can be 
constrained surprisingly well. The bias uncertainty does 
not significantly demolish the promise of this aspect. In 
the Wq-Wi panel, the two EOS parameters wq and W\ are 
strongly correlated. This suggests that the redshift range 
0.2-0.4 is not large enough to distinguish evolved model 
and non-evolved model, and that the constraints on w 
mainly depend on an effective value of w. In fact, the de- 
generacy is represented by wq + 0.13wi = const., which 
means the effective EOS is given by w e s = w(z = 0.13). 
Since the EOS affects the apparent clustering through 
Eq.(|l]), it is reasonable that the effective redshift is 
halfway to the sample volume. 

To see how the uncertainty of the EOS evolution af- 
fects the constraints on f^Qo, w c r, we re-parameterize the 
EOS by w(z) = w c ff + (z — 0.13)wi. In Fig. ||, expected 
bounds are plotted for the case when one of the dark- 
energy parameters Oqo, w e fi, u>i ar e marginalized over. 
The error bounds on Qq and w e g are barely affected by 
w\ marginalization. This means that our method is not 
sensitive to the uncertainty of the EOS evolution, and is 
only sensitive to the effective value of the EOS. 

In the above analyses, the cosmological parameters 
other than that of the dark-energy component and bias 
parameter are fixed. These other parameters can be de- 
termined from various cosmological observations. In or- 
der to estimate the effects of the uncertainties from these 
other parameters, error bounds between dark-energy pa- 
rameters and other parameters are plotted in Fig. 0, 
without any other marginalization. Each panel in this 
figure indicates a correlation between each pair of pa- 
rameters. Correlations with erg ar e not plotted, but are 
quite similar to those with bias 6, as naturally expected. 
The parameter fiQo is positively correlated with J7ko an d 
h. This is because the galaxy clustering accurately con- 
strains fl-Moh = (1 + SIko — ^Qo)h. If we do not use 
prior knowledge of either the curvature or the Hubble's 
constant, the error bound on Oqo is enlarged by a fac- 
tor of 4-5. There are negative correlations between each 
dark-energy parameter and the bias parameter. This can 
be understood by noticing that both parameters increase 
the amplitude of fluctuations when the present amplitude 
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FIG. 2: Expected error bounds of dark-energy parameters 
with the effective EOS w e g — w(z — 0.13). One of three 
parameters are marginalized over. 
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The prior input of the curvature and the Hubble's con- 
stant are important to reduce the bounds on the dark 
energy density, but not so important on the EOS param- 
eters. 
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FIG. 3: Error bounds between dark-energy parameters and 
other parameters. Solid line: wo = — 1 model. Dotted line: 
wo — —0.5 model. 



is fixed. The EOS parameters are quite independent on 
the curvature and Hubble's constant. 

In summary, we predict for the first time expected 
degree of constraints on dark energy models from 
intcrmediate-rcdshift galaxies such as LRGs, and show 
that £Iqq and w e s can be constrained surprisingly well. 
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